Pair tunneling through single molecules 
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By a polaronic energy shift, the effective charging energy of molecules can become negative, 
favoring ground states with even numbers of electrons. Here, we show that charge transport through 
such molecules near ground-state degeneracies is dominated by tunneling of electron pairs which 
coexists with (featureless) single-electron cotunneling. Due to the restricted phase space for pair 
tunneling, the current-voltage characteristics exhibits striking differences from the conventional 
Coulomb blockade. In asymmetric junctions, pair tunneling can be used for gate-controlled current 
rectification and switching. 
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Introduction. — Electronic transport through single 
molecules is distinguished from transport through quan- 
tum dots by the coupling to well-defined vibrational de- 
grees of freedom. This coupling has two principal con- 
sequences: First, it leads to the emergence of vibra- 
tional sidebands in the current-voltage {IV) character- 
istics. This phenomenon occurs at large voltages which 
exceed the vibrational frequency, oj. Its theoretical de- 
scription was pioneered in Refs. 0,0 and extended in, 
e.g., Refs. nil 13, following experiments [10,1101 on 
transport through molecular junctions. 

Second, the coupling to molecular vibrations induces a 
polaron shift and can lead to a negative effective charging 
energy U. In physics, the concept of negative-C/ centers 
was first pointed out more than three decades ago [lo| . 
and is realized in many amorphous semiconductors. In 
chemistry, the scenario of negative U is known as "poten- 
tial inversion" [1^ . An important ingredient in realizing 
this scenario in molecular junctions may be a reduction 
of the true molecular charging energy by screening due 
to metallic electrodes or an electrolytic environment 
[l3| . Unlike vibrational sidebands, negative U manifests 
itself in transport through molecules already at low volt- 
ages. The prime manifestation of negative U studied to 
date [ij, list II6I llTl llSj concerns the Kondo transport at 
very low temperatures. 

In general, finite on-site interaction (of any sign) opens 
a new transport channel between molecule and leads. 
Namely, at finite U, two electrons can hop onto the 
molecule simultaneously. But it is only for negative 
U that this process can dominate the transport. In- 
deed, negative U favors even electron occupation num- 
bers of the molecule. A Kondo resonance occurs when 
two states, whose occupation numbers differ by two, are 
degenerate. Then, a virtual pair transition assumes the 
role of the spin- flip leading to the formation of 

Kondo correlations in the ground state. However, the 
fact that the relevant pseudospin degree of freedom is 
associated with charge makes this Kondo state rather 
fragile. The underlying reason for this was elucidated by 
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FIG. 1: (a) Level configuration of the negative-f/ model. 
Shown are the one-particle energies for the singly-occupied 
dot (sd), the doubly-occupied dot (sd + U/2), and the energy 
for which holes can propagate through the doubly-occupied 
dot {sd + U). The right panel illustrates the relevant types of 
processes: (b) cotunneling, (c) and (d) pair tunneling. 



Haldane , who pointed out that for the charge Kondo 
effect, any deviation from degeneracy acts as a Zeeman 
field on the pseudospin and, thus, suppresses the Kondo 
correlations. This observation has a drastic consequence 
for transport. For the conventional (spin) Kondo effect, 
the low-temperature conductance eventually reaches the 
unitary value over a wide range of gate voltages. By 
contrast, for the charge Kondo effect the unitary value 
is only achieved precisely at the resonance gate voltage. 
This fragility of the charge Kondo effect in transport 
through molecules with negative U was illustrated by 
several recent numerical simulations [Hill [13. The 
main message of the present paper is that negative U 
strongly affects the transport through a molecule even 
in the rate-equation regime at high temperatures, where 
Kondo correlations are irrelevant. We show that in this 
regime, a negative U leads to a unique scenario for the 
passage of current through the molecule, which we study 
analytically at all gate voltages and biases. 

Effective Hamiltonian. — Consider a molecule with a 
single spin-degenerate single-particle orbital Sd dominat- 
ing the transport between two metallic electrodes. As 
is sweeped by a gate voltage, the electron occupation Ud 
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of the molecule switches from empty (n^ = 0) to doubly 
occupied {ud = 2) at the resonant gate voltage defined 
by the condition 2ed + U = (We measure all en- 
ergies from the Fermi energy.) Due to the negative U, 
single occupation (rid = 1) of the molecule is unfavorable 
at any gate voltage. A schematic of the corresponding 
configuration of one-particle energies near resonance is 
shown in Fig. ^a). Sequential tunneling of single elec- 
trons is clearly exponentially suppressed at bias voltages 
\eV\ <C \U\. Instead, the dominant sequential transport 
processes near resonance involve the coherent transfer of 
electron pairs Representative pair-tunneling pro- 

cesses are illustrated in Figs, ^c) and (d). It is impor- 
tant to realize that in addition to processes in which the 
two electrons enter the molecule from the same lead, the 
electron pair on the molecule can be created (or anni- 
hilated) by two electrons tunneling in from (or out to) 
opposite leads. In parallel with these sequential pair- 
tunneling processes, single-particle cotunncling processes 
[see Fig. Qlb)] also take place. We find that while both 
processes are of the same order, they can be easily distin- 
guished since the single-particle cotunncling contribution 
does not exhibit any structure near the resonance gate 
voltage. 

To capture all eight pair-tunneling processes sys- 
tematically, we perform a Schrieffer- Wolff transforma- 
tion We start from the Hamiltonian H = 
Hraoi + -ffvib + -ff loads + ^^i Q, H , whcrc the sepa- 
rate contributions correspond to the electronic molec- 
ular orbital i/moi = + Uudiridi, to the molec- 
ular vibrations H^ih = hwh% + Xfnoihi^ -)~ b)nd, to 
the leads iJicads = Ea=i,fl Ep,^ ep4p^Cap^> and to 
the tunneling between leads and molecule. Hi = 
J2a=L,RJ2p,a {taclp^da + h.c.) . Here, the operator da 
annihilates an electron with spin a on the molecule, Capa 
annihilates an electron in lead a {a = L,R) with mo- 
mentum p and spin ct, vibrational excitations are an- 
nihilated by h, and t^ ji denotes the tunneling matrix 
elements. The electron-phonon coupling (with coupling 
constant A) can be eliminated as usual by a canonical 
transformation 0, 0| , which implies renormalizations of 
the tunneling Hamiltonian ta tae^''^^^^ of the or- 
bital energy Ed — A^fiw, and of the charging energy 
U ^ U — 2X^fiuj. It is this last renormalization which 
opens the possibility of a negative effective charging en- 
ergy. The underlying reason for this renormalization is 
that the energetic polaron shift is proportional to the ex- 
cess charge of the molecule squared which is of the same 
form as the Coulomb charging energy. To be specific, we 
focus on temperatures and biases where only the vibra- 
tional ground state is populated. In this case, we obtain 
an effective Hamiltonian iJcff = ^^moi + -ffioads + Hi with 
negative effective charging energy U and Franck-Condon 
suppressed tunneling matrix elements ta taS^^ 1"^ (ab- 
sorbed into the definition of ta in the following). 

Performing the Schrieffer- Wolff transformation on _ffoff 



in the usual way |22| , we eliminate Hi to lowest order and 
obtain a transformed Hamiltonian 



-ffsw — -ffmol + H\ 
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where we ignore a term which merely renormalizes Sd- 
For positive U, one only retains 
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~t~ tata' ^^i^^ap) ^^(T^CT^apcT^a'p'cr C^apa^o-' 1?' (^^^da^ h.C. 

which describes the direct and exchange interactions be- 
tween molecule and lead. Here, we have introduced the 
abbreviations M(e) — - jz^i^^zfuj^ and f, I =i,T- 
With U > and in the absence of charge fluctuations, 
transport predominantly proceeds via single-electron co- 
tunneling described by -ffdir.ex- By contrast, for negative 
U it is crucial to retain the pair-tunneling terms 

iJpair = tata' M{eap)d^diclp,^clp^ + h.C. (3) 

aa'pp' 

Obviously, i/pair contributes only for nonzero effective 
charging energy U . 

Rates. — Using Fermi's Golden Rule, the pair Hamilto- 
nian -ffpair leads to the rate 

WS^, = I de M\e)fa{e)fa' {2ed + U^e), (4) 

for pairs of electrons tunneling onto the molecule. Here, 
Ta = ini'a \ta\^ IS the energy scale of single-particle tun- 
neling in junction a for a constant density of states Va- 
The superscripts {aa') denote the leads from where the 
spin-up electron (a) and the spin-down electron (a') orig- 
inate. The analogous rates 14^2^0 ^'^^ pairs leaving the 
molecule are obtained from Eq. by replacing each 
lead Fermi function fa with a factor (1 — fa). 

In the regime where single-particle occupation of 
the molecule is negligible, i.e. \2ed + U\ ,\eV\ , ksT ^ 
£d, \£d + U\, the integral in Eq. Q approximately re- 
duces to an integral over the Fermi functions alone, 
and can be expressed in terms of the function F(e) = 
e/[exp(/3e) — 1]. For symmetric voltage splitting [i^, i.e., 
/a(e) = /(e — eVa) with Vl^r = ±V^/2, the explicit result 
for pair tunneling reads 



^^M^{0)F[±{2ed + U + eVa + eVa')], (5) 



2 (x = 



where the upper (lower) sign refers to x = 
2 ^ 0). 

The pair-tunncling rates, Eq. jsj, have several remark- 
able features, (i) Pair tunneling with electrons originat- 
ing from the same lead [a = a', Fig. GJc)] leads to bias- 
dependent rates. By contrast, pair tunneling with elec- 
trons from different leads [a ^ a'. Fig. ^d)] is bias in- 
dependent, as the energy missing in the lead with lower 
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Fermi energy is exactly compensated for by the addi- 
tional energy available in the lead with higher Fermi en- 
ergy, (ii) When pair tunneling is energetically allowed, 
its rate is proportional to the detuning of the pair state 
2ed + U from the relevant Fermi energy. This is in 
sharp contrast with rates for single-electron sequential 
tunneling, which are independent of the energy of the 
single-particle state. This unusual behavior of the pair- 
tunneling rates arises from the fact that only the sum 
of the energies of the two tunneling electrons is fixed 
by energy conservation, making the phase space for pair 
tunneling proportional to the detuning [2^. We will see 
below that this leads to characteristic features of pair 
tunneling both in the gate-voltage and bias dependence 
of the conductance. 

Similarly, the rates for cotunneling from lead a 
to a' can be obtained from -H"dir,ox as W'gi^.Q = 
2TaTa'F{eVa)/he'^, including an explicit factor of 2 for 
spin. The corresponding rates W2^2 are obtained by the 
replacement (ed + U)"^. 

Having established the relevant processes and their 
rates, we now describe transport outside the Kondo 
regime by the corresponding rate equations [3 Issf . 
Since the occupation probability for the singly-occupied 
molecule is negligible for \eV\ \U\, the stationary rate 
equations reduce to = i-2W^2^o — PqWo^2, with the so- 
lution Po = W2^o/[W2^o + Wo^2l P2 - 1 - Pq. Here, 
Wi^f denotes the total rate for transitions from initial 
state i to final state /, i.e., Wx — J2a a' l^x" • The sta- 
tionary current / = Jp^"^ + 1"^°^^ evaluated in the left 
junction, involves contributions from pair tunneling and 
cotunneling, 



p-"ye = PoWo^2-P2W2^o, 

r°'/e^ P0V0^0+P2V2^2- 

Here, the coefficients v and w are given by 



(6) 
(7) 
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(8) 

The factor of 2 in the last equation originates from the 
coherent transfer of two electrons in this pair-tunneling 
process. [All remaining coefficients are obtained from 
Eq. (jH)) by interchanging "2 <-> 0" in the subindices.] 

Results. — The rate equations provide a complete ana- 
lytical description of the nonlinear current-voltage char- 
acteristics through negative-[/ molecules. Specifically, we 
find for the linear conductance 
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+f{-2ed-U)^ + f{2ed + U) 



Here, the last two terms in brackets, which arise from 
cotunneling, give a slowly varying background. By con- 



11; 



1 ' 1 ' 1 
(a) 4 

ii 


- 


1' ' 

~ t\\ 








— T = 


5.0 X 10--' \Ul 


--■ T = 


2.5 X 10"^ |c;i 


T = 




1.0 X 10-= \U\ 
1 , 




-C.04 -0.02 0.02 0.04 

ed. + U/2 {\U\) 



— eV = 
-- eV = 3 X 10-^ \U\ 
eF = 6 X 10-= \U\ 

I I I I I I I I I 

-0.04 -0.02 0.02 0.04 

ed + U/2 {\U\) 



FIG. 2: (color online) Conductance as a function of gate volt- 
age, based on our analytical results, (a) Linear conductance 
G for several temperatures. The conductance curve exhibits a 
featureless background due to cotunneling and a distinct peak 
of constant height and width ~ T due to pair tunnehng. (b) 
Conductance for several bias voltages at T = 1.0 x 10""^ \U\ 
for a symmetric junction. With bias voltages eV S> fcsT, the 
width of the double peak is given by eV . 



trast, pair tunneling, described by the first term in brack- 
ets, leads to a remarkable conductance peak. As illus- 
trated in Fig. Ufa), its height Gq = 2Ae'^V lT R/U'^h ex- 
ceeds twice the cotunneling background and is temper- 
ature independent, while its width is proportional to T. 
This feature, which is a direct consequence of the phase 
space for pair tunneling, should provide an unambigu- 
ous experimental signature of pair tunneling as opposed 
to ordinary Coulomb blockade conductance peaks, whose 
integral strength is temperature independent. 

While the linear conductance is identical for symmetric 
and asymmetric devices, their finite-bias behavior differs 
considerably. For large voltages \eV\ » fcsT, we find for 
the pair tunneling current 



^pairs _ 



2TlTr 



2TLrR \s\ 



(10) 



e+TL\S+\{TR\S^\+TL\5\) + {L^R^ 



where S = 2ed + U , 5± = (5 ± eV , and 0± = 
0{-5)9{'^5±) - 9{5)0{±5±). For symmetric devices, Eq. 
()10|l implies that the width of the conductance curve is 
fixed by the bias voltage, as illustrated in Fig.|2tb). This 
behavior is in stark contrast to the conventional Coulomb 
blockade, where one finds two symmetric sharp peaks at 
detunings ±V/2 from resonance. This difference which 
also arises from the phase space for pair tunneling, is 
further emphasized in Fig.Ofa) for all bias and gate volt- 
ages. 

Devices with a significant asymmetry in the molecule- 
lead couplings exhibit a striking asymmetry with respect 
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FIG. 3: (color online) Current as a function of bias and gate 
voltage for (a) a symmetric junction with Tl ~ — ksT, 
(b) an asymmetric junction with Tl = O.lfcsT, Fr — lOksT, 
and with ksT = 0.5 ■ 10""^ \U\. For asymmetric junctions, one 
observes the rectification effect due to pair tunneling. 

to voltage inversion, see Fig. |3Ib). We find in this case 
that pair tunnehng causes strong current rectification, 
whose transmission direction can be switched by a gate 
voltage. To understand this rectification effect, consider 
an asymmetric device with Tfj ^ Tl, and suppose that 
the Fermi level is higher in the left lead. Then, the pair- 
tunneling current proceeds as follows: (i) A pair with 
electrons from opposite leads jumps onto the molecule; 
(ii) the electron pair on the molecule is transferred to the 
right lead. While (ii) is a fast process (~ F^), the cur- 
rent is limited by process (i) with rate ~ T^Tji. Then, 
the switching by a gate voltage immediately follows from 
the rate for process (i), which is exponentially suppressed 
only for 2ed + U > 0. By a similar analysis for the oppo- 
site bias, one finds an exponential suppression of the pair 
current for 2ed + U < 0, thus also explaining the recti- 
fication. These features render molecules with negative 
U interesting candidates for devices with transistor-like 
characteristics. 

Conclusions. — Near the degeneracy point, transport 
through a molecular junction with negative U differs 
drastically from the conventional Coulomb-blockade sce- 
nario at positive U. For negative U, current flow makes 
the molecule alternate in time between even occupation 
numbers, which is accomplished by pair tunneling. This 
is the dominant transport mode in the gate-voltage - 
bias- voltage domain shown in Fig. ^ 

It is intriguing that the on-site attraction of two elec- 
trons makes pair tunneling through molecules qualita- 
tively similar to tunneling through a superconducting 
grain, considered in Ref. |26]. However, in the case 
of molecules, the physical picture of transport is more 
complex. This is because a pair can be created on the 
molecule by electrons tunneling from different leads, as 
illustrated in Fig.^d). By contrast, for a grain with size 
larger than the superconducting correlation length, elec- 
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FIG. 4: Stability diagram for the ground-state occupation nd 
and dominant transport modes as a function of gate and bias 
voltages for a symmetric device with negative U. 

tron pairs enter and exit the grain only from the same 
lead, i.e. only the processes in Fig. ^c) are responsible 
for the passage of current [8(ij |. As an important con- 
sequence of this difference, negative- J7 molecules act as 
gate-controlled rectifiers. 
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